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Abstract 

This paper extends the quantum search class of 
algorithms to the multiple solution case. It is shown that, 
like the basic search algorithm, these too can be repre- 
sented as a rotation in an appropriately defined two 
dimensional vector space. This yields new applications - 
an algorithm is presented that can create an arbitrarily 
specified quantum superposition on a space of size N in 

0(jN) steps. By making a measurement on this super- 
position, it is possible to obtain a sample according to an 
arbitrarily specified classical probability distribution in 

O(jN) steps. A classical algorithm would need £l(N) 
steps. 
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1. Introduction 

(i) Sampling: Sampling is a fundamental technique for 
approximating answers that cannot be directly, or effi- 
ciently, computed. The uses of sampling for gathering 
information is the primary domain of the field of statis- 
tics. Sampling has also been linked to complexity the- 
ory, e.g counting, volume estimation, primality testing. 
For some applications uniform sampling is enough but 
for most applications one needs sampling according to a 
specified probability distribution. It is known how to 
rapidly generate samples when the probability distribu- 
tion has a certain structure, e.g. when the distribution is 
log concave, samples can be generated in logarithmic 
time [1]. In general, the probability distributions are 
only specified implicitly, i.e. it is possible to calculate 
the probability of any given point easily but there is no 
known structure to the distribution. It is not possible to 
sample precisely according to an arbitrary probability 
distribution in fewer than Q.(N) classical steps. For 
example, it is easily seen that we will need to examine at 
least half the points, for if we leave out half of the points 
we could be missing a point with an arbitrarily high 
probability, e.g. the situation where all the probability is 
localized in a single point. 

A quantum mechanical system is in multiple 
states simultaneously, quantum computing algorithms 
such as quantum search make use of this feature. This 
paper extends quantum search to develop a quantum 
algorithm that can create an arbitrarily specified super- 



position (and hence an arbitrary probability distribution) 
on N points in only 0(*jN) steps. 

(ii) Quantum computing: Just as classical digital sys- 
tems can be constructed out of two state systems called 
bits, quantum mechanical systems can be constructed 
out of basic two state quantum mechanical systems 
called qubits. Quantum mechanical operations that can 
be carried out in a controlled way are unitary operations 
that act on a small number of qubits in each step. A good 
starting point to think of quantum mechanical algo- 
rithms is probabilistic algorithms [2] (e.g. simulated 
annealing). In a quantum mechanical algorithm, the sys- 
tem is started in a state that is easy to prepare on which a 
sequence of simple operations is applied. When the sys- 
tem is observed after applying these operations, it gives 
the answer to a difficult computational problem with a 
high probability. 

Just like classical probabilistic algorithms, quan- 
tum mechanical algorithms work with a probability dis- 
tribution over various states. However, unlike classical 
systems, the probability vector does not completely 
describe the system. In order to completely describe the 
system we need the amplitude in each state which is a 
complex number - such a specification is called an 
amplitude vector or a superposition. The probabilities in 
any state are given by the square of the absolute values 
of the amplitude in that state. The evolution of the sys- 
tem is obtained by premultiplying this amplitude vector 
by a state transition matrix, the entries of which are 
complex in general. It can be shown that in order to con- 
serve probabilities, the state transition matrix has to be 
unitary, i.e. the columns of the transformation matrix are 
orthonormal [2], It is easily shown that the following 
three transformations are unitary: 

(i) NOT - a one-input one-output gate. The output is 



the inversion of the input. 

(ii) CNTRL-NOT - a two-input two-output gate. The 
first output is the same as the first input. If the first input 
is 1, the second output is the inversion of the second 
input; if the first input is 0, the second output is equal to 
the second input. 

(iii) CNTRL-CNTRL-NOT - a three-input three-output 
gate. The first two outputs are equal to the first two 
inputs respectively. If the first two inputs are both Is, 
the third output is the inversion of the third input; if 
either of the first two inputs is 0, the third output is equal 
to the third input. 

These are the analogs of classical NAND and NOR 
gates. Using these three gates it is possible to synthesize 

any boolean function f(x) that can be synthesized clas- 
sically with approximately the same number of gates. In 
order to develop more powerful quantum mechanical 
algorithms, we need some operations that are basically 
quantum mechanical, i.e. the entries in the state transi- 
tion matrix are not all 0's and l's. The two such opera- 
tions that we will need in quantum search are the Walsh 
Hadamard (W-H) transformation and the selective 
inversion of the phase of certain states, these are dis- 
cussed in the following two paragraphs. 

A basic operation in quantum computing is the 
operation M performed on a single qubit - this is repre- 



sented by the following matrix: M=-^— 

J2 



1 1 
1 -1 



, i.e. 



the state is transformed into a superposition where the 

two states & 1 have the same amplitude of , this 

J2 

superposition is denoted by: -^=(|0) + |1)) ; similarly 
state 1 is transformed into the superposition 



— (|0)-|1)). A system consisting of n qubits has 

72 



N=2 basis states. If we perform the transformation 
M on each qubit independently in sequence, the state 
transition matrix representing this operation is of dimen- 
sion NxN. Consider a case when the starting state is 
one of the N basis states, i.e. a state described by a gen- 
eral string of n binary digits composed of some Os and 
some Is. The result of performing the transformation M 
on each qubit will be a superposition of states consisting 
of all possible n bit binary strings with amplitude of 

n 

~2 1 

each state being ±2 , i.e. ±— — . This transformation is 

Jn 

referred to as the W-H transformation [8]. 

The other transformation that we need is the 
selective phase inversion of the amplitude in certain 
states. The transformation matrix describing this for a 4 
state system with selective phase inversion of the sec- 



ond state is: 



10 
0-100 
10 
1 



Unlike the W-H transforma- 



tion, the probability in each state stays the same. A 
realization of this kind of transformation can be 
achieved using the gates discussed so far [4]. 

In the exhaustive search problem, a function 

f(x) defined over N states (denoted by x) is given 

which is known to be non-zero at a single value of x , 

say t (t for target) - the goal is to find t . If there was no 

other information about f(x) and one were using a clas- 

N 

sical computer, then on the average it would take — 

function evaluations to solve this problem successfully. 
[9] found a quantum mechanical algorithm that took 



only 0(jN) steps. 

The quantum search algorithm [9] consisted of 

Jn repetitions of the operator -I.WI-W starting with 

the state (here W denotes the W-H transformation, I ( 
denotes the selective phase inversion of the target state 
t, I- denotes the selective phase inversion of the 

state). Later on it was discovered that similar results are 
obtained by replacing the W-H transform by almost any 

quantum mechanical operation (say U) and the state 
by any basis state s. It was shown that by starting with 

1 



the basis state s, and carrying out O 



repetitions 



of the operation sequence -IJJ I f U , one could reach 

the t state [10] (similar results are also proved in [5]). 
This showed that one could use any starting state and 
unitary operation U and from these amplify the ampli- 
tude in a desired target state t. A new class of algorithms 
was thus invented. These extended far beyond search 
problems - in fact, it was shown that this framework 
could be used to enhance almost any quantum mechani- 
cal algorithm [10]. 

One constraint with these algorithms was that 
they worked when the problem has exactly one / state. 
For many problems like game-tree search, this was a 
major restriction and the algorithms either could not be 
shown to work with multiple solutions or needed com- 
plicated workarounds [6]. The paper [11] mentioned the 
multi-solution case as an open problem. 

2. This paper This paper shows how the quantum 
search algorithm with general unitary transformations 
can be extended to the multi-solution case. From this, an 



algorithm for generating an arbitrarily specified super- 
position is developed. This can be used to sample 
according to a general probability distribution - the 
number of steps required is approximately the square- 
root of that of the corresponding classical algorithm. 
The following is the organization of the rest of this 
paper: 

section 3 carries out an analysis of the generalized 
search algorithm with an arbitrary number of solu- 
tions, 

section 4 considers the case with a single solution, 
section 5 extends this to the multi-solution case, 
section 6 shows how this can be used to generate an 
arbitrarily specified superposition. 

3. Generalized search (arbitrary 

number Of Solutions) The algorithm starts 
with the system in an s state (s for source), the object is 
to drive the system into the t states, note that there is a 
single s state but multiple t states. The following sec- 
tions show how this can be accomplished by means of 
an operation Q which is defined as the composite opera- 
tion Q = -I s U i I { U (note that, following standard 
matrix notation, this denotes the following sequence of 
operations: first U, then / , then U 1 and finally -/ .) 
/ is a diagonal matrix with all diagonal elements 

equal to 1 except the ss' h elements which is -1 . In 
Dirac notation, which is commonly used for matrix 
operations in quantum mechanics, / may be written as 

I s = I-2\s)(s\. Here / is the identity matrix and \s) 

the column vector with all except the s' h element equal 

th 

to zero, the s element is 1 , (s\ denotes the corre- 



sponding row vector. Similarly I f is a diagonal matrix 

th 

with all diagonal elements equal to 1 except the tt 
elements which are equal to -1 . I f can be written in the 

form I t = I- ^2|f)(f| . U denotes an arbitrary unitary 

t 

matrix and U 1 is its inverse. The following results 
hold for arbitrary U , section 6 describes how to choose 
U based on the specified probability distribution. 
In this notation: 

Q\s) = -(I s U l I t U)\s) 



= -(I-2\s){s\)U^ I~Y}\t){A 

= -\s) + 2\s) + ^{t\U\s)U l \t) 

t 

-^\s){s\u-\t){AU\s) 



U\s) 



Note that (t\U\s) = U u . Also, since U is unitary, U is 

equal to the transpose of the complex conjugate of U , 

-1 * 

and therefore (s\U \t) = (t\U \s) which is equal to 
U ts . Therefore the above equation becomes: 

q\s) = \s)+Y J ^u ts )u Y>- y, a \ u a 2 \ s) 

t t 

Similarly, for any t state: 
QU~ l \t) = -(I s U l I t U)U l \t) 



= {I-2\s){s\)U^\t) 



trV>- 2u u \ s ) 



The analysis so far shows that, if there are r| t states, 
then the (r| + 1) dimensional space defined by the vec- 
tor \s) and the r| vectors U \t) is preserved by the 
operator Q . The next two sections show that there is 
indeed a simpler two dimensional subspace that is also 
preserved by Q . 

4. Single solution Case First of all consider 
the situation with a single t state. This is the same as the 
generalized search situation considered previously [10]. 
In this section, we quickly rederive the main result 
since, in the next section, this will help to analyze the 
multi-solution case. The following two equations follow 
from the previous section: 

Q\s) = \s)(\-A\U t ^) + 2U ts U' l \t) 



QU V> = U l \t)- 2U* \s) 



Thus Q is a transformation in the 2-dimensional vector 

space defined by the two vectors \s) and U 1 |f) ■ ln 
matrix notation this may be written as the following 
transformation: 







JJ-\t) 





{l-A\U ts \ ) -2U ts 
2U ts 1 



In order to find the effect of repeated applications 
of Q , we use standard matrix analysis which consists of 
finding the eigenvalues and eigenvectors of the transfor- 
mation matrix. Assuming U f to be real and small, the 
two eigenvalues and eigenvectors are approximately: 



X, = l+2iU ts , Vl = 



& 



X 2 = 1- 2iU ts , v 2 = 



The initial state vector is |s) , which in terms of 
the eigenvectors may be written as i(vj + v 2 ) . After r| 
applications of Q , this transforms into: 
^(vjA-j + VjXj) which may be simplified to: 



COS(2f/ fJ Tl) 

sin(2f/ ;s ri) 



Therefore, 



41/.. 



applications of Q, 



transform |s) into U \t) . From this, \t) may be 
obtained by a single application of U . After this, a mea- 
surement will reveal the target state, f, with certainty. 

5. Multi-SOlution Case The analysis of section 
4 breaks down if there is more than one t state with dif- 
ferent values of t/ This section extends the previous 
results. The relevant equations from section 3 are: 



Q\s) = \s) 



QU V> = U V>- 2U*\s) (for each t state.) 



Multiply the second equation by t/ and sum over all t 
states. The equations now become - 



Q\s) = \s) 



1 _ 4 £ I U tsf + 2 Y, U t* U 1 |f) (aS bef ° re) 



-1 



-1 



Q is thus a transformation in the two dimensional com- 



L-l 
u ts u \t) . 

t 

L-l 
U fs U \t) and denoting 



l^l^l by u , the above transformation can be rep- 



resented as: 



Q 



\s) 



(1-4k ) -2m 
2u 1 



This is exactly the same transformation as in the 
previous section with a single t state. Therefore, as in 



the previous section, if u « 1 , then — applications of 



2 transform the vector \s) into -\ U ts U \t) . The 

f 

case when u is (9(1) is similar but not particularly 
interesting since this can be solved with a high probabil- 
ity by just applying U to \s) once and carrying out an 
observation. 

Note that the algorithm in the previous section 
was only shown to work when U f was real. This sec- 
tion shows that a very similar result holds in the single 
solution case, even when 17 is complex. 

6. State vector engineering Suppose that we 

are required to synthesize a specified superposition. The 
amplitude in each of N states, denoted by x , is required 
to be proportional to a given complex valued function 
f(x) . Since f(x) can be scaled by an arbitrary factor, we 



assume, without loss of generality that the maximum 

value of \f(x)\ is less than or equal to 1 . 

The following algorithm synthesizes the speci- 



fied superposition in approximately ^ — steps. 

4 f- — • i - ,2 



One immediate application of this algorithm is in sam- 
pling according to an arbitrary probability distribution 

in O(jN) steps. For this a quantum superposition is 
generated and then through a measurement in an appro- 
priate basis, a sample is obtained. 

Solution Define N=2 n states by n qubits, denoted 
by x . Include an additional ancilla qubit. Initialize the 
state so that all qubits are in the state - the state of the 
whole system is denoted by: (0, 0) (the first denotes 
a single qubit in the state and denotes each of the 
other n qubits in the state.) 

Next consider the following unitary operations 
(which constitute the building blocks for our algorithm): 
t/j : leave the first qubit unaltered and apply a W-H 

transform to the other n qubits. 
U 2 ■ carry out a conditional rotation of the first qubit so 

that the state: (0, x) gets transformed into the fol- 
lowing superposition: 

(f(x)(0, x) + Jl -\f(x)\ 2 (l, x)) , i.e. the amplitude 
of the state (0, x) is f(x) and the amplitude of the 



state (l,x) is Jl - \f(x)\ . This type of unitary 
operation has previously been used in quantum 
computing algorithms, e.g. in the mean estimation 
algorithm in [10]. It can be accomplished by first 



transferring f(x) into the phase through conditional 
phase inversion and then converting it into ampli- 
tude information. 
I f : in case the first qubit is , invert the phase; if the 

first qubit is 1 , leave it unchanged. In other words, 
states with the first qubit in the state are t states. 
/ : in case all the qubits (including the first qubit) are 

0, invert the phase; else leave it unchanged, i.e. 
(0, 0) is the s state. 
Clearly if we start with the (0, 0) state and apply 
U l and then U 2 , the amplitude in the (0, x) state is 

m 



. In other words, if we define the composite opera- 



tion U= U 2 U l , the s state as (0, 0) and the t states as 
the (0, x) states; then {/ , the matrix element between 

/(*) 



s and the relevant t state is 



It immediately follows from section 5, that by 
starting with the (0, 0) state, and applying the sequence 



-1 , 



of operations defined by Q = -U S U I ( U), 



n N 



4 l£l/wl 2 



times, followed by a single application of 



U , we get the first qubit in the state and the remain- 
ing n qubits in a superposition with the amplitude of the 

fix) 



x state as 



£[f(*)f 



(i) Let f(x) be 1 at T) points in the domain and zero 

everywhere else. Since ^|/"(x)| 2 = r| , the algo- 

x 

% In 



rithm needs - /— steps. After this it reaches the 
4A/r| f 

same superposition as reached in the basic quan- 
tum search algorithm with r| solutions and it 
needs exactly the same number of steps as the 
quantum search algorithm to reach this. 

(ii) The number of steps required by the algorithm 

2 

depends on ^ \f(x)\ . In case this quantity is not 

X 

known in advance, the superposition can still be 

synthesized in O(jN) steps, by trying out the 
algorithm with a few carefully chosen runtimes. 
After this, the ancilla qubit is measured - the 
algorithm is repeated until this is observed to be 
(it can be shown that, with appropriately cho- 
sen runtimes, the probability of not getting even a 
single falls exponentially with the square of 
the number of times the procedure is repeated 
[7]). Once a is observed, the remaining n 
qubits immediately collapse into the desired 
superposition. 

(iii) The algorithm assumed that the function f(x) 
was normalized so that its maximum value was 
equal to 1 . The algorithm as presented in this 

section is equally valid for different f(x) , pro- 
vided the maximum value does not exceed 1 . For 
example, if the desired probability at each x is 

— 2 _ ^ 

specified (i.e. \f(x)\ ), the value of ^[f(x)| 



7. Observations 



becomes 1 and the algorithm needs exactly 
^ JN iterations to attain this distribution. In gen- 
eral, the number of steps required is smaller if we 
can choose a larger constant to scale f(x) . This is 
maximized for the choice made in this section 
where the maximum value of f(x) is 1 . 

(iv) The analysis presented above is very similar for 
the situation where the initial state \s) , instead of 
being a basis state is an arbitrary superposition as 
considered in [3]. The main difference is that in 
the computational basis, I s (the matrix that 

inverts the phase of the source state) will not be 
diagonal. Of course, it will be diagonal in any 
basis that includes \s) . In terms of the Dirac nota- 
tion it will still be described as I $ = I-2\s){s\ , 

i.e. it reflects the \s) state and leaves all orthogo- 
nal states unchanged. 

(v) A superposition with multiple target states, of the 
type discussed in section 5, is a pure state in a dif- 
ferent basis. Thus the algorithm rotates from a 
given source state to another pure state. This is 
just like quantum search except that the final state 
is a basis state in a different basis - the quantum 
search algorithm could only rotate into a pure 
basis state in the same basis. 



algorithm. In addition to the basic result that a quantum 

computer can search a domain of size N in 0(jN) 
steps, its contribution has been to inspire several new 
ideas and algorithms. This paper presents the most 
recent such development which shows that the quantum 
search class of algorithms can be made to work in the 
presence of multiple target states. Such a result may 
seem obvious for classical computers; however, for 
quantum computers one needs to take into account inter- 
ference effects due to the different solutions and the sit- 
uation is more complicated. 

An application has been presented where one can 
generate a sample according to an arbitrary probability 
distribution in a number of steps which is only a square- 
root of that required by a classical algorithm. Another 
immediate application is in extending the framework of 
the quantum search class of algorithms with arbitrary 
unitary transformations of the type discussed in [5] & 
[10], so that they are no longer limited to problems with 
a single solution. 

Just like previous versions of the quantum search 
algorithm, the framework here is completely general 
and it is expected that it will find more applications. 
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8. Conclusion The quantum search algorithm is 
perhaps the simplest possible quantum mechanical algo- 
rithm that yields a significant advantage over a classical 
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